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I. 



INTRODUCTION 



In writing Problem Oriented Computer Languages several 
non numerical applications of computer programming must be 
mastered and the literature about such activities is nearly 
nil for the PORTRAN language. It is to these problems that 
this work addresses itself. The professional programmer 
would probably elect to write these codes in machine 
language but the user who only knows FORTRAN cannot afford 
time to learn the various machine languages in use today. 
All of the features discussed below are programmable in 
FORTRAN aad illustrations are given in each use. 

The CAL language is then -used to illustrate how these 
features are incorporated in a problem oriented language. 
The CAL program was developed by Professor Edward L. Wilson 
of the University of California, Berkeley, California, for 
Structural Mechanics and Structural Engineering. 

CAL combines matrix manipulation routines with direct 
stiffness computation options to produce a program for the 
automated analysis of structures. CAL also has significant 
capabilities as an instructional tool in linear algebra. 

In order to use any structural analysis program it is 
necessary to idealize the structure into a series of joints 
connected by structural members. The joints are commonly 
referred to as nodal points and the structural members as 
elements. The program input consists of: 

a) nodal point locations; 

b) information on how nodes are connected; 

c) boundary conditions that are applied to the 

nodes ; 



d) properties of the elements; 

e) the loading to be analyzed. 

The program output, at a minimum, consists of nodal point 
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deflection values and normally includes nodal forces and/or 
element stcess values. CAL breaks this procedure down into 
a series of simple steps under the control of the user. It 
was designed to execute rapidly on small problems, so 
students can quickly see results of an analysis. There are 
options which allow the user to debug data without printing 
previously obtained results. Looping operations are also 
available. This allows a user to program iterative 
numerical algorithms, greatly expanding the usefulness of 
the program. In this version, CAL will solve problems with 
about 50 degrees of freedom (DOF) . 

The primary aim of the ‘thesis work presented here was 
to modify CAL for use on the IBM 360/67 computer at the 
Naval Postgraduate School (NPS) . The NPS version of CAL has 
both interactive and batch operating modes, and is nearly 
machine independent. 

Most of this thesis effort was not involved in the 
numerical procedures ased for ca Iculations, but rather in 
the non-nuaer ical areas of data input and processing. For 
this reason, numerical procedures will not be discussed. 
Subsequent chapters will provide: 

a) a substantive discussion of some non-nuraerical 
techniques useful in writing scientific programs such as 
CAL ; 

b) a discussion of the internal organization of 

CAL-NPS; 

c) instructions for the implementation of CAL on a 
different computer. 

The appendices contain a listing of the FORTRAN code and an 
Instruction Manual for CAL-NPS. 

Modifications to the program were carfully made to 
avoid changes in the instruction manual written by Professor 
Wilson. PeQ commands have been added and three commands 
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have been extended for more flexibility in the interactive 
use of the program. The original instruction manual has 
been modified to reflect these changes and is reproduced in 
Appendix 3. The author gratefully acknowledges Professor 
Wilson's permission to use the manual. 

The non- numerical techniques discussed in Chapter II 
are this author's effort to explain how several 
state-of-the-art programming techniques work. It is hoped 
that readers will find this section beneficial in their own 
programming efforts. The CAL organization discussed in 
Chapter III is an excellent outline for general scientific 
programs. Additionaly, Chapter III provides details for 
writing user supplied subroutines. For readers who are 
interestel only in how to use CAL, Chapter II and IV may be 
omitted . " . 
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II. 



NON-NUMERICAL TECHNIQOES IN PORTRAN 



Chapter II is used to discuss FORTRAN programming 
techniques for writing general scientific programs. The 
reader is assumed to have some familiarity with the FORTRAN 
language and computer programming (i.e. , completed CS2700 
FORTRAN Programming or equivalent at NPS) . After 
introductory remarks on the attributes of a scientific 
program, dynamic dimensioning is discussed in Section A, 
free and object time formats are discussed in Section B, 
data management is discussed in Section C, and finally, 
looping is discussed in Section D. The terms will be 
defined as encountered in these sections. 

A general purpose program for the computer solution of 
a class of scientific problems should be simple to use, 
flexible and reliable. This requires a program that: 

a) is modular, i.e., sections of program may be 
deleted, inserted or modified with ease; 

b) has numerical methods that are stable and accurate 
for the type of problem encountered; 

c) may require the user to have a knowledge of 
numerical methods necessary for the selection of the 
appropriate solution technique or boundary conditions but 
does not require knowledge of computer programming beyond 
the ability to prepare the input data for a program, 
initiate execution of a program, and create data files; 

d) is not limited by array sizes in dimension or common 
statements ; 

e) permits iteration, i.e., the repeated execution of a 
user selected sequence of operations; 

f) requires a minimum of central processor time and 
internal computer memory space (core) for execution; 
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A very useful method for meeting these requirements is 
to structure the program around the use of a problem 
oriented language. Such a language consists of a set of 
mnemonics representing each possible step in the general 
solution of the problem. 

A mnemonic is a word that serves as a memory aid to 
identify some thing or some action. In FORTRAN, the 
mnemonic ZDS is used to represent the process of computing 
the cosine of an angle. All mnemonics can be thought of as 
representing a function subprogram, which can be very 
complex. The use of mnemonics to represent various 

operations in the formulation and solution of a particular 
problem is a very attractive method of writing a program for 
the general solution of a class of scientific problems. 
This technique promotes modularity and greatly simplifies 
use of the program. The mnemonic used for a particular 
operation or computation should be representative of the 
operation. For instance, a programmer might select the 

mnemonic SDL VS for the operation of computing the vector x 
from the matrix equation Ax = B where A and B are known. 
The mnemonic identifying a particular operation together 
with a list of arguments necessary for the computation is 
read from one data card. This card is followed by a series 
of data cards to provide data, not currently in core 
storage, for computation. The ability to enter the 
information in free format or to allow object time format 
specification dramatically reduces the work required to use 
the program. 

References (1), (2), and (3) provide examples of 

problem oriented language programming. Apart from these few 
examples, there is very little written on the non-numerical 
techniques needed to write a program for the general 
solution of a scientific problem using a problem oriented 
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demonstra: es this process graphically. Usually it is 
possible to use a compiler that will optimize the machine 
language cole to minimize the run time for a program in 
execution. Optimizing compilers take more time and require 
more core space to compile than do non-optimizing compilers. 
It is possible to compile and link-edit a program then store 
the load module on a mass storage device (magnetic tape or 
disk) . This always results in saving central processor 
(CPC) time for subsequent program runs and can result in a 
significant reduction of the amount of core required. 



'FORTRAN 

f Program 
/statements/ 




/PROGRAM; 
Oi/rpirr / 



Program Processing 
FIGURE 1 
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The alternative to the classical style of programming 
is dynamic dimensioning. In order to understand the concept 
of dynamic dimensioning, it is necessary to have some 
appreciation for how a computer handles a common statement. 
Inside the computer, the common statement will cause the 
reservation of sufficient bytes of memory to store listed 
variables. A byte is a segment of a computer word which can 
contain either one alpha (hollerith) character or about two 
decimal digits. If the same common requires a different 
number of bytes in different subroutines, the largest 
requirement will be the one reserved. Any type of variable 
maf occupy the reserved storage as long as the address of 
the starting byte of the variable is consistent with the 
type of variable. Different types. of variables frequently 
have a different number of bytes in a word. For instance, a 
complex variable will normally have twice as many bytes per 
word as a real variable. To illustrate this point, consider 
a program where all variable names beginning with A are 
double precision (8 bytes) , all variable names beginning 
with Z are complex double precision (16 bytes), and 
variables beginning with N are integers (4 bytes). The named 
common EXAMPLE might be listed as follows in different 
subrout ines : 



In each case the same number of bytes of storage are 
reserved (48) but the common has a significantly different 
type of variable! Introductory FORTRAN courses normally 
teach the student that a particular common should always be 
an exact duplicate each time it appears in the program 



COMMON /EXAMPLE/ N ( 1 2) 
COMMON /EXAMPLE/ A (3, 2) 
COMMON /EXAMPLE/ Z(3) 



( 1 ) 

( 2 ) 

(3) 
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without »uch other discussion of the uses of common 
statements. The example lines above contrast sharply with 
cautious style but can prove very useful in actual 

I gives the address of the starting 



this 

problem solving. Table 



byte for each variable in EXAMPLE. 



BYTE 

1 

5 

9 

13 

17 

21 

25 

29 

33 

37 

41 

45 



VARIABLES 

,A (1,1) ,Z (1) 
,A(2,1) 

,A (3,1) , Z (2) 



7 , A ( 1 , 2) 

9[' ,A (2,2) , Z (3) 

1 J) 

\l ' M3 ' 2 > 



Byte Addresses for Variables 
TABLE 1 



Note that N(3), A (2,1) and the imaginary portion of Z (1) all 
have the stime byte as their starting address. Additionally, 
48 bytes would have been reserved even if the array 
dimension of both (2) and (3) had been set at one since (1) 
still requires 43 bytes. One more feature of importance is 
that A ( 2, 1 ) , A (1,2) and A(3,2) contain the same information 
as the coresponding imaginary portions of the Z vector. This 
fact is quite useful in programming with complex variables. 

The common EXAMPLE above demonstrates how a computer 
stores a two dimensional array but amplification may prove 
helpful. The right most subscript in a multiply dimensioned 
array is always the one that increases slowest as the array 
is going into computer storage. The address in memory of 
an element of a linear array is computed much faster than 
the address of an element from an array with multiple 
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subscripts. The address of the element A(I,J) of a two 
dimensional array can be computed in a corresponding linear 
array as follows: 

L = ( ( J- 1) *N R ♦ I - 1 ) * I PR + 1 (4) 

where 

SR is the number of rows in the matrix 
IPR is the ratio of bytes per word of the matrix to 
bytes per word of the linear array. 

For example, compute the index in N of A (3,2) from the 
common EXAHPLE: 

here 1=3, NR=3, J=2 and IPR = 8/4 = 2 
L = ( (2-1) *3 ♦ 3 - ‘ 1) *2 * 1 = 11 
Compare this answer with the results in table 1. 

Dynamic dimensioning is a process that utilizes the 
computer's ability to store different types of variables in 
the same common array. Array storage is reserved in blank 
common with a large, one-dimensional, integer array 
(hereafter called the main array) for all subscripted 
variables needed for a particular computation step. This is 
normally done in a vary short main program. Problem size can 
be controlled by the dimension of the main array in the main 
program. The dimension of the main array in all subroutines 
is set a; 1 and never changed. Typically, to change a 
problem size, two cards of a six card main program need to 
be altered. See figure 2 for an example. The main program 
needs to be compiled, but there is no need to use an 
optimizing compiler. The remainder of the program can be in 



15 



a load module on a mass storage device. 



C MAIN PROGS AM 

COMMON MAX,L(6000) 

MAX = 6000 
CALL DRIVER 
STOP 
END 

C SUBPROGRAM TO CONTROL EXECUTION 

SUBROUTINE DRIVER 
COMMON MAX, L ( 1 ) 



CALL COMPUTE (L (N1) , L (N 2) , . . . , N3 , N4) 



RETURN 

END 

C SUBROUTINE FOR COMPUTATION 

SUBROUTINE CO MP UT E (A , B , . . . , NROH,ICOL) 
DIMENSION A (NROW, 1) ,B (1) 



RETURN 

END 

Dynamic Dimensioning Sample Program 



FIGURE 2 



The main program calls a subprogram which controls 
execution, reserves storage in the main array for 
subscripted variables, and provides the address of array 
arguments in the main array to the computation subroutines. 
This subprogram must also ensure that there is sufficient 
room in the main array before a subscripted variable is 
allowed to be stored. Boundary misalignment occurs if the 
type of variable being stored requires more than one integer 
wocd for storage, and the starting address of the 
subscripted variable in the main array is not an address 
that is compatible with allowable addresses for the variable 
being stored. Some compiliers automatically correct for 
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boundary alignment but often this is accomplished with a 
considerable penalty in the time required for execution. 
Consider, for example, that a real double precision variable 
on an ,I3M 360 computer requires two integer words for 
storage. A boundary misalignment will occur unless the 
double precision variable is at an odd integer word address 
counting from the first integer word in blank common. In 
figure 2, MAX is the first integer word and L(1) is the 
second word in blank common. L(2| is, therefore, the first 
address available for the storage of the double precision 
element. L(2) and L(3) are required for the storage of one 
double precision number. If the double precision variable 
started iu L ( 1 ) , it would be out of proper alignment for 
this type of variable. A simple way to compute a proper 
starting iddress is: " • 

SI = NSIZE + MOD ( N VAR * NSIZE , IPR ) + 1 

where 

N 1 is the starting address to be computed. 

SVAR is the number of integer words used by the 
variables in blank common preceeding the main array, L. 

NSIZE is the number of words in the main array 
that have been used. 

IPR is the number of integer words required to 
store the variable being processed. 

MOD (I,J) is an ANSI standard function which 
returns the remainder from ( I / J ) . 

A method to keep track of the starting addresses for 
subscripted variables will be discussed later. The amount 
of storage available may be compared with the amount 
required for storage of a subscripted variable as follows: 
LEFT = MAX - ( N1 + NEL * IPP. ) + 1 

where 

LEFT is the number of integer words remaining. 
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array. 



MAX is the number of integer words in the main 



N1 is the starting address in the main array for 
the variable. 

NEL is the number of elements in the variable. 
Note that when blanks are left between the end of an array 
and the start of the next array to compensate for boundary 
alignment, NEL of the leading array is effectively changed. 

IPR is the number of integer words required per 

el ement . 

Positive steps must be taken to prevent LEFT from 
becoming negative. One method to compute the number of 
words that have been used in the main array is 

N SI ZE = BAX - LEFT 

where all variables are used as defined above. 

Figure 2 is an example of some programming necessary to 
implement the dynamic dimensioning concept. The variable 
MAX and the dimension of the L array must have the same 
value in the main program. Outside the main program, the L 
array neads only to be dimensioned one since adequate 
storage is already reserved. In the control subprogram, N1 
and N2 are used to provide starting addresses for 
subscripted variables to the computation subroutine. N3 and 
N4 are used to provide array dimensions. The dimension 
statement in the computation subroutines can use variables 
for array subscripts. Linear arrays and the final subscript 
of multiply subscripted arrays need only be set to one as 
shown in figure 2. 
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B. FREE AND OBJECT TIME FORMATS 



Recall that in classical programming, data entered into 
the computer must be entered exactly in accordance with the 
format statement that has been compiled with the program. 
For instaice, suppose that the coordinates and a reference 
number are required for a point in 3-space in 
format (3F1 5.5,11 0) . Figure 3 is the data card required for 
the point ( 1 . 0,- 1 .0, 1 . 0) assigned reference number 3. Note 
that if tie 3 is not in card column 55, the reference number 
assigned to the point will be wrong. When the input device 
is a terminal, the probability of making an error on 
entering data with this type of fixed format is high. 
Errors in entering data frequently require the user to 
restart tne program and re-enter all the data to correct a 
small mistake. 



1 



FIGURE 3 

Object time format (also called variable format) 
permits the user to enter the format for subsequent input or 
output operations as a part of the data deck at the time the 
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program eaters execution. This means that it is possible to 
read an array prepared in (F10.0) format and an array 
prepared in (5E15.7) format without changing the format 
statement and recompiling the program. 

Object time formatting is accomplished by reading a 
card containing the format for the data to follow, then 
reading or writing the data as specified in the format, A 
sample of the required programming to read data is: 

DIMENSION FOR (20), M3), ••• 

• 

1 00 READ (5,200) FOR 
200 FORMAT (20A4) 

300 READ (5, FOR) ( A (I) ,1=1 , 3) , J 



Statement 100 represents an ordinary read command in FORTRAN 
according to the format statement numbered 200. This format 
statement causes 80 columns of hollerith data to be read (4 
columns ei ch into 20 variables). Statement 300 causes 4 real 
numbers to be read into the vector A. However, instead of 
being read according to a numbered format statement 
contained in the program, these numbers are read according 
to the information contained in the variable FOR. FOR might 
contain : 

(3F4.0,I1) 

causing the coordinates of the point discussed above to be 
read as the vector A in format (3F4.0) and the reference 
number to be read as the integer J in format (II). Figure 4 
shows the cards required to enter the point and reference 
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Sample Object Time Dimension Data 



FIGURE 4 



The first non-blank character on the data card containing 
the format must be an open parenthesis. The last non-blank 
character must be a close parenthesis. The characters in 
between must constitute a valid format for the data to be 
read or written. 

Object time dimensioning allows the user to adjust 
formats to suit the data being input. However, if the user 
neglects to enter a format or enters one containing an 
error, the resulting problems can be frustrating. Free 
format relieves the user of the difficulties associated with 
using any type of predefined format. Figure 5 is an example 
of the data card required to enter the point and reference 
number discussed earlier. This example uses a comma (,) as 
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an end of field symbol and a blank ( ) as the end of record 
sya bol . 
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Format Data 



FIGOBE 5 

The first step of programming free-format input is to 
read the data card into an array called the card-image 
array. The data card must be read in format ( 8 0 A 1 ) so that 
the Jth element of the card-image array represents the 
character in column J of the data card. It is important to 
note that the value stored in the computer for a digit read 
as a hollerith character is different from the value stored 
for the same digit read as a number. The card-image array 
is then processed element by element to form the input 
variables. These input variables may be integers, real 
numbers or character strings, however, the variables must 
still be in the order expected by the program. If the 
expected input is a data card containing a step name and 
matrix argument list, a card containing real numbers will 
cause an error. Free format input allows the program and 
not the system error handling routine to determine 
appropriate corrective action. 
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The first step in processing an element of the 
card-imagj array is to identify the character stored in the 
element. 3dm precise the identification must be depends on 
the application. The character being used as an end of 
field symbol and the character being used as an end of 
record oust alMays be identified precisely. If, for 
instance, the expected data is a step name and matrix 
argument list, then the programmer might chose to identify 
all alphabetic characters as one type of symbol, all digits 
as a second type of symbol and all characters not otherwise 
identified as a third type of symbol. However, if the 
expected data contains a real number that may or may not 
have an exponent the choice of characters might be; 0, 1, 2, 
3, 4, 5, 5, 7, 3, 9, * , ., D, E, and others. The variable 
NSIM will be used to indicate the type of character. 

The history of the processing '.is important in deciding 
what to do with a character. For instance, in the 
assembling of the real number 1.21, there are two types of 
character, digits and a decimal point. If the character 
encountered is a digit, the proper assembly action depends 
on whether the decimal point has occurred. If the decimal 
point has not occurred, the numbec being assembled must be 
multiplied by 10 and then the digit added to the number. On 
the other hand, if the decimal point has occurred, the digit 
must be divided by 10 then added to the number being 
assembled. The power of 10, N, is a counter representing 
the number of places to the right of the decimal this digit 
belongs. The varialbe IACT will be used to indicate the 
appropriate assembly action. 

In addition to integral and fractional digits, a number 
may carry an algebraic sign and/or a sign in the exponent. 
The sign lay occur at the beginning of the number or at the 
start of the exponent, however, a sign encountered elsewhere 
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will defiaetly be an error. It follows that for the 
assembly process of real numbers, the program performs one 
of the following operations: (a> sign identification, (b) 

integral digit assembly, (c) decimal point identification, 
(d) fractional digit assembly, and { e ) exponent assembly. 
These processes are refered to as states. If the state is 
integral ligit assembly and the character encountered is a 
sign, then obviously an error has occurred. On the other 
hand, if the character encountered is a decimal point, then 
there is io error but the state for the next operation needs 
to be changed to fractional digit assembly. & variable named 
ISTATE is introduced to indicate which one of the above 
operations the program last performed. Now, continuing with 
the assemoly of a real number, the variable ISTATE may be 
assigned as follows: 



ISTATE Processing: 

1 at the beginning of a Field. 

2 after the beginning but before the 
decimal point. 

3 after the decimal point. 

4 after the exponent sign. 

5 after an error has occurred in the 
field. 

Again, for different applications, the programmer may define 
other states as needed. 

The next step in processing a character from the 
card-image array is two-fold. The proper action to assemble 
the variable must be executed. This, of course, is given by 
the value of IACT. Also, if required, the state must be 
changed. A single integer can be used to represent both the 
next state and the assembly action required. For example, 
the inforiation contained in ISTATE and IACT can be combined 
into a single variable, IPROC, by: 
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IPEOC = I ACT* 10 0 + ISTATE 
The inforiation can be separated by: 

I ACT = IPHOC/100 

and 

IS TATE = MOD (IPROC, 1 00) 

A matrix can be assembled that has the row index defined by 
1ST ATE and the column index defined by NSYM . If an element 
of this matrix represents IPROC for the ISTATE and NSYM 
that define the element, then this matrix is an efficient 
method of programming the assembly task. The matrix 
constructed is called the Symbol-State Matrix. Refere 
(4) contains futher discussion of Symbol-State matrices 
their uses. 

For illustration, consider the problem of reading 
real linear array of length NEL. For simplicity, eleme 
of the array may not have exponents. The state 
identified as: 

ISTATE Processing: 
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d) Identify the character occupying the ICARD column of 
the card-image array (NSYM) from the hollerith translation 
table. 

e) Using NSYM and the current ISTATE as indices for 
entry in the Symbol-State matrix, called ISST, obtain IPROC 
(IPR3C = ISST (ISTATE, NSYM) . From IPROC, compute IACT and 
the new ISTATE. 

f) Carry out the appropriate assembly action based on 
the value of IACT as follows: 

1) Set a flag indicating the number is negative. 

2) Assemble the digit iato the integral part of 
the number (X = 10.0 * X + DIGIT). 

3) Assemble the digit into the fractional part of 
the number (X = X * DIGIT/ (10. 0**NCOUNT) ) , and step the 
fractional digit counter (NCOUNTj . 

4) Print an error message. An action to 
compensate for the error may be included here. 

5) Store the element in the storage array. If the 
array is filled, return; otherwise, reset NCOUNT and 
continue. If the symbol causing this action is EOR, reset 
ICARD ana continue at step B. 

Note that a plus sign or decimal point changes the 
state of the program but does not reguire an assembly 
action. After the assembly action above is complete, proceed 
to step 6. 

6) Increment ICARD and test for the logical end of 
record. If ICARD is greater than the record length, reset 
ICARD, execute step f.5 and continue at step b. If ICARD is 
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l e3S than or equal to the record length, continue at step d 
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on the value of IACT. 
if there is no decimal 
assumed after the last 
Symbol-State Table or 
difficult once the cone 
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string, for instance, 
bytes per integer word 
array may contain: 



This table does not indicate an error 
point. Instead, the decimal point is 
digit of the field. Additions to the 
modifications to the logic are not 
ept is understood. 

problem of assembling a character 
as a matrix name. If there are U 
, then one element of the card-image 
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vhere "bl" represents the blank character. The byte 
containing the hollerith character S from the card-image 
array must be moved into the appropriate byte of the 
character string variable being assembled. This is not a 
trivial problem. ANSI FORTRAN does not support byte 
manipulation. Major computer manufacturers have added 
features to their FORTRAN which allow byte manipulation for 
string assembly. CDC allows for byte manipulation with BASK 
and SHIFT operations. In IBM FORTRAN the job can be done 
with a sunroutine as follows: 



SUBROUTINE SHI FT ( A , B , I , J) 

LOGICAL* 1 A(1), Bill 

THIS SUBROUTINE MOVES THE JTH BITE OF 
VARIABLE B INTO THE ITH BITE OF 
VARIABLE A. 

A (I) = B (J) 

RETURN 

END 



The next update of the ANSI FORTRAN Standard tentatively 
includes variable word size specification with a 
"CHARACTEl *n" declaration, where n is an integer. If the 
change materializes then the above subroutine will work for 
any machine by changing LOGICAL*! to CHARACTER*!. 
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C . DATA MANAGEMENT 



The comments on data management have been restricted to 
the management of data storage in core. Only one of several 
alternatives is discussed and the management system 
presented is not necessarily the most efficient method. It 
is advocated because it represents a straight forward 
approach and one does not have to be a computer scientist to 
implement the procedure. 

The problem of in core data management can be reduced 
to the problem of storing, finding and deleting arrays in 
the main array mentioned previously. One way of organizing 
the main array is to store an- array directory at the start 
of each array in the main array. The array directory needs 
to contain the following information: 

a) The number of elements in the array. 
b| The precision of the array (integer 
wor ds/elenent ) . 

c) A unique identification for the array. 

d) The number of rows and/or the number of columns 
in the array. 



D1 Ml 



D2 



M2 



D3 



M3 



Diagram of Array Storage 

The process of reserving storage for an array should: 

a) Ensure adequate space is available in the main 
array. Processing must be interrupted for resolution of the 
problem if sufficient room is not available. 

b) Adjust NSIZE to compensate for the additional 
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array . 



c) Create and store the array directory based on 
information from the calling program. 

d) Check boundary alignment and adjust NEL for the 
previous array if necessary. It is possible to adjust NEL of 
the current array so that any subsequent array will oe on a 
proper boundary no matter vhat the precision of the 
subsequent array. This is a simple solution but does make 
some available space unavailable. 

e) Return the main array address for the first 
element of the array to the calling program. 

Locating an array requires searching the array 
directories for a match with the identifier of the array. 
One method of doing this type of matching is: 

IF (COMP (L (N) ,11) ) 30 TO nnn 

where 

N is the address, in the main array, of the first 
word of the array identifier being tested. 

LI contains the name to be matched. 

ann is the statement number that is branched to if 
a match is found. 

COMP is a function subprogram defined as follows: 

FUNCTION COMP (I , J) 

DIMENSION I (3) , J (3) 

LOGICAL COMP 
COMP = .FALSE. 

DO 10 K = 1,3 

IF (I (K) . NE. J (K) ) RETURN 

. 10 CONTINUE 

COMP = .TRUE. 

RETURN 

END 

This example assumes that three integer words are used to 
contain the array identifier. If the array is found, the 
main array address of the starting element and the directory 
informatic i needs to be returned to the calling program. If 
the array is not found, processing should be interrupted for 
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Skyline storage schemes make more efficient use of 
storage by simply not storing many of the matrix elements 
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that are zero. The skyline storage algorithm is the Host 
storage efficient of the two types. Beference (1) contains 
a substantive discussion of the skyline storage Method, and 
examples of matrix manipulation subroutines for use with the 
skyline algorithm. Suffice it here to say that the skyline 
storage alogrithm requires the creation of two vectors. 
One is a vector of the elements from the matrix that are 
actually stored. The second vector contains the indices in 
the first vector of the diagonal elements of the matrix. 
For this case, it is convenient to use the same directory 
for both vectors. If there are NEL elements under the 
skyline, each element takes IPR integer words and there are 
NR rows in the matrix, then the storage required is 
NDIR * NR + NEL * IPR 

a typical block of the main array might look like 



DIRECTORY DIAGONAL ELEMENT INDICES 



STORED ELEMENTS 



The above presents a straightforward approach to data 
management in core. It is codeable in ANSI FORTRAN IV by an 
individual familiar with the basics of the language. It 
does not require a professional programmer for 
i mplementation. 

D. LOOPING 

The capability to repetitively execute a series of 
instructions (looping) is an important tool in problem 
solution. Frequently, for non-linear problems, there is no 
direct method of solution. An iterative scheme is normally 
employed to solve this type of problem. Brute force 
iteration is possible with any program by making multiple 
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runs of the program with suitable alteratio 
data for each run. This is very wasteful of 
time and computer time. On the other hand, 
possible, the user can program an iteration s 
commands available, then place the series of 
loop. & relatively simple method of provi 
capability is to store the commands of the 1 
in the main array. Then, with appropriate 
possible to execute a series of instructio 
from memory, just as if they were coming 
device. Details of the method of loopin 
program CiL are contained in’Chaptar III. 

A particularly attractive advantage of 
ability to iterate for selected eigenvalues 
This abiLity permits one to 'approximate 
significant digits in the computer soluton o 
follows : 

NSD = NDA - ALOGIO (EMAX/EMIN) 

where 
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NSD is the namber of significant digits. 

NDA is the number of digits available in the word 

size. 

EMAX is the maximum eigenvalue of the matrix. 

EMIN is the minimum eigenvalue. 

ALOGIO is the ANSI standard function for the 
logarithm to the base 10. 

This approximation holds for positive, definite, symmetric 
matrices. Computer results for problems where NSD 
approaches 0 should be viewed with extreme skepticism. 
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III. CAL ORGANIZATION 



This chapter provides a general overview of the 
operation and organization of CAL. It is intended to 
provide sufficient information to permit users to program 
subroutines for use in CAL and to operate the CAL system 
seasibly. 

CAL is a modern computer program written in the FORTRAN 
I V programming language. It was written as a teaching aid 
for illustrating the direct stiffness method of structural 
analysis. In addition to extensive linear algebra 

capabilities, the program provides several analytical 
alternatives for both static and dynamic analysis of elastic 
structures. 

CAL execution is flexible, controlled by user selection 
of operations in a logical sequence (not unique) from the 
operations available in the program. The input data deck is 
a sequence of modules. The first card of each module 

contains the name of the requested operation, a list of 
matrix argument names, and integer parameters used in that 
operation. Comments may be placed on the card following a 
blank. 

The operation card looks like this: 



OPERATION, Ml, M2,. . . , M9 , N 1 , . . . N4 COMMENTS 
Ml througi M9 are matrix argument names. A matrix argument 
name consists of up to 8 alpha-numeric characters, the first 
being alphabetical. N 1 throuch S4 are positive integers. 
An operation card may use from zero through nine matrix 
arguments and/or from zero through four integers. If 
required, the operation card is followed by a sequence of 
data cards containing numerical information used for the 
step. Some features of the output from CAL are also 



34 



controlled by user selectable operations. Details of the 
operations available and how to execute CAL at NPS are 
presented in Appendix B. 

CAL is logically divided into four segments. The main 
segment contains the main program and service subroutines 
used commonally by all other segments. In subsequent 
descriptions, the main segment will be referred to as CAL. 
Among other things CAL contains the subroutines to deliver 
the next operation and the subroutines for data management. 
The three remaining segments will be referred to as GROUP 1, 
GRDUP 2, and GROUP 3. 

GROUP 1 contains subroutines which perforin matrix 
manipulations, i.e., load, print, multiply, etc. GROUP 2 
contains subroutines associated with formulation of the 
static problem and display of fhe computation results. For 
example, there are subroutines to: • 

a) input nodal geometry; 

b) input boundary conditions,' 

c) input loading conditions,* 

d) form element stifness and mass matrices by 
several different methods; 

e) combine element matrices into global matrices; 

f) display results. 

Finally, GROUP 3 contains subroutines associated with 
dynamic analysis. Two methods of evaluating the second 
order egiations of motion, eigen analysis, and printer 
plotting routines are available. 

Matrix storage in CAL is done dynamically as described 
in Chapter II. It is not possible to state with precision 
the maximum problem size. Experience shows that with the 
default parameters of the NPS version, static analysis 
problems are limited to approximately 48 degrees of freedom 
(DDF). Dynamic analysis problem size is controlled by both 
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the system DOF and the number of 3 ispl acemen t vectors to be 
calculated. For a 25 DOF system, about 100 displacement 
vectors can be calculated. 

All arrays and matrices in CAL are stored in a one 
dimensional matrix specified in blank common as L(XXXXX). 
XXXXX represents the total storage reserved for subscripted 
variables. Each array is preceeded by a directory which 
contains: 

1) the number of elements in the array (HEL) . 

2) the number of integer words reguired to store 
an element (IPR) , 

3) the number of columns in the array, 

4) the number of rows in the array, 

5) the first four characters of the matrix name, 
and 5) the last four characters of the matrix name. 

An array uses (NEL*IPR *■ 6) integer words of storage. 
Larger problems can be solved by increasing the size of the 
L array in the main program. The variable MAX must be set 
to the current dimension of the L array. This data 
management system is described in Chapter II. CAL uses 
three subroutines for the data management sub-system. There 
is a subroutine called LIST which is used to reserve storage 
for new matrices as they are created. The calling arguments 
are, in order, NM, NR, NC, IPR, and IERR. NM is an array 
name assigned by the user. Array names must start with an 
alphabetic character and conatin 1 to 8 alpha-numeric 
characters. 

Since the name, Nfi, is used only in the directory, 
there is no reason to conform to FORTRAN name conventions. 
Thus, names beginning with the characters I through N may be 
assigned for real matrices with impunity. Call the 
stiffness matrix K, the mass matrix BASS and the nodal 
coordinate array NODES. All variables are treated as real. 
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double precision numbers unless clearly specified as 
integers in the command description. 

NR and NC are respectively, the number of rows and the 
number of columns for the matrix to be stored. IPR is the 
number of integer words required to store an element. IEBR 
is an error parameter. If IEBR is not equal to 1, then an 
error has occurred during the list operation. 

The subroutine LOCATE is used to locate arrays in the L 
array. The calling arguments are, in order, NM , NA , NR, NC, 
and IERR. NM, NR, NC, and IERR are used as described above. 
NA is the index in the L array for the first element of the 
array. Ii a subroutine call, L (NA) is used to pass the 
matrix NM to and from the called subroutine. 

The remaining data management subroutine is DELETE. It 
has a single calling argument. NM is used as discussed 
previously. DELETE removes the matrix NM from the L array 
and frees the storage taken by NM at the end of the L array. 
Figure 7 Illustrates this process. 
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Now consider once more the function of the main program 
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segment CkL. In addition to the data management subroutines 
described above, the CAL segment contains four subroutines 
used to deliver the next operation to the calling program 
segment, and a timing subroutine. The subroutine INPUT 
determines whether the next operation is part of a loop or a 
new command from the input device. If the next operation is 
not part of a loop, INPUT calls the subroutine OPREAD. In 
turn, OPREAD calls the subroutine RCARD to read and 
interpret the free format information of the next command. 
Pree format data input is discussed in Chapter II. If the 
command initiates looping, then OPREAD calls RCARD until all 
operations in the loop have been read. As operations are 
received by OPREAD, this subroutine stores them row-wise in 
the L array. After the last operation of the loop has been 
read, OPREAD calls LIST to prepare the directory. INPUT 
calls OPSTOR to provide each command in sequence from the 
loop matrix in the L array until' the loop operations are 
complete. The subroutine SECOND is called by INPUT to 
determine the amount of time used in each computation step. 
A block diagram of the connections is: 




A black diagram of the total system follows: 
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UP 1 contains two operation names 
routines. These names are USERA 
subroutine might be used, for 
uivalent nodal forces for a 

eeded' for a user operation. The 
is coded, normally in a subroutine 
sed in the CAL program. A buffer 
itten to interface the numerical 
data management system. To 
interface, consider the familiar 
ation (already available in CAL) . 

to compute C where C = A*B. The 
multiplication is 



C(I,J) = A (I. K) *B (K,J) 

A subroutine for the job might be: 

SUBROUTINE MULTIPLY (A, B, C, NRA ,NCA, NCB) 
DIMENSION A (N R A , N CA) ,B(NCA,NCB) ,C(NRA,NCB) 
DO 20 I =1 , N RA 
DO 20 J = 1 , N CA 
X = 0.0 
DO 10 K=1 , NC A 
10 X = X + A (I,K) * B (K, J) 

20 0 (I,J) = X 
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RETURN 

END 

The opera: ion card night be selected as: 

USERA, Ml r M2, !13 COMMENTS 

where HI, M2, and M3 are respectively the matrix argument 
names for the dummy arguments A, B, and C. The recommended 
steps to be programmed in the buffer subroutine named USERA 
are : 

a) locate Ml and M2; 

b) ensure the number of columns in A is equal to 
the number of rows in B; 

4 

c) call the DELETE subroutine to ensure that a 
matrix with the same name as M3 does not already exist in 
storage ; 

d) call LIST to reserve storage for M3; 

e) test TERR to ensure that no errors have 
occurred in the above steps; 

and f) call subroutine MULTIPLY to complete the 
operation ; 

At the minimum, steps a, d, and f must be programmed. Here 
is a possible subroutine to meet the requirement. 

SUBROUTINE USERA (IERR) 

COMMON MAX ,NDP,L (1) 

COMMON /CARD/ IN HOL (3 , 1 0 ) , N (4 ) 

CALL LOCATE (INHOL (1,2) , IA, NR A ,NCA,I ERR) 

CALL LOCATE (INHOL (1,3) , IB , NRB ,NCB, I ERR) 

IF ( NR A . NE . NCA ) IERR= 2 
CALL DELETE (INHOL (1,4) ) 

CALL LIST (INHOL (1 ,4) ,NR A , NCB , NDP , IE RR) 

CALL LOCATE (INHOL (1,4) , IC , NRA , NCB, I ERR) 

IF (IERR .NE. 1 ) RETURN 

CALL MULTIPLY (L(IA) , L (I B) ,L (IC) , NR A , NCA ,N CB) 

RETURN 
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END 



LOCATE, LC3T, and DELETE are used as described previously. 
Blank common appears so that OSERA has access to the L array 
and the parameter ND?. The variable N DP is initialized at 
2, the number of integer words required to store a double 
precision variable. The named common CARD is initialized by 
HCARD or OPSTOR. The INHOL array contains the operation 
name and the matrix argument names. Each name is stored 
columnwise, i.e., column 1 has the operation name (USERA in 
this case*, column 2 has the name of thematrix argument Ml, 
etc. Row 1 elements have the first four characters, and row 
2 elements contain the remaining four characters of the 
respective argument. Row 3 is not used in the NPS version 
of CAL. In the first call to subroutine LOCATE, INHOL(1,2) 
provides the name assigned by the user to matrix argument 
Ml. On return from LOCATE, IA will be the index in the L 
array of the first element of Ml. IB and IC respectively 
provide the same information for M2 and M3. Other variables 
are used as discussed previously. 

Now, generalizing the procedure of the above example, 
the user must program two subroutines to define a user 
operation for CAL. One subroutine performs the operation, 
and the second acts as a buffer between CAL and the 
operation. Two names are available, OSERA and OSERB, for 
the buffer subroutine. Both have a single passing argument, 
the parameter IERR. The program assumes an error has 
occurred when IERR takes on a value other than 1. The 

buffer suDroutine must: 
a) 



operation, 



b| 



locate matrix arguments designated for 
reserve storage for arrays created by 



the 



the 



operation, 



c) pass the elements of the L array for the matrix 
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parameters used in the operation. 

In addition, it is recommended that the buffer subroutine: 

a) ensure that the matrix arguments to be used are 
compatible with the operation, i.e., square, symmetrical, 
etc . 

b) call delete with names of matrices to be 
created prior to reserving storage. This will prevent 
matrices with duplicate names. 

c) test TE8R to ensure no error has occurred prior 
to calling the operation subroutine. 

The operacion subroutine can be coded with normal PORTRAN 
t echn iques. 

In this chapter, the internal organization of CAL has 
been described, along with several important features of the 
program. Most of the non-numerical procedures described in 
the preceeding chapter are illustrated in CAL. This chapter 
has highlighted the flexibility of CAL and provided an 
example to assist users in coding their own operations to be 
executed by CAL. 
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IV. 



IMPLEMENTING CAL ON A DIFFERENT COMPUTER 



In this chapter, the procedure to convert CAL for use 
on a different computer system is discussed. The NPS 
version of CAL has been organized to minimize the difficulty 
associated with implementing the program on a new computer. 
There renins, however, some things that are machine 
dep endenr. 

The following information is required for the new 
system . 

a) How many characters can be stored in one 
integer word? 

b) How many integer words does it take to store a 
real variable? For structural, analysis, a real variable 
should contain at least 12 significant figures. 

c 

c) What method can be used to manipulate bytes? 

d) How is elapsed CPU time computed? 

The answers to the above questions determine the complexity 
of the conversion task. 

A. IBM FORTRAN STATEMENTS 



There are a few statements in this version of CAL 
solely in the IBM FORTRAN Language. The statement "IMPLICIT 
REAL*8 (A-H, 0-Z) " makes all real variables double precesion. 
This statement appears in almost every subroutine. It 
should be replaced by an appropriate statement when using 
another computer. The calls to ERRSET suppress printout of 
error messages. These statements appear together in 
subroutine CALI and must be removed. They do not need to be 
replaced. The call to SETIME in the main program should be 
removed for use on non-IBM machines. 
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B. TIKE COMPUTATION 



There is no standard method of computing elapsed CPO 
tine. The subroutine SECOND must be rewritten to compute 
elapsed C?J time in seconds for the new system. The library 
of subroutines for the computer system should have a 
suitable routine. 

C. BYTE MANIPULATION 



As previously mentioned there is no standard way to 
shift bytes in a variable. The subroutine KOVEB will work 
if one byte logical variables ar-e available. If this is not 
the case, than MOVEB must be .rewritten using a byte 
manipulation technique suitable for the computer system. 

D. OPERATION AND MATRIX NAMES 

In this version, names are stored four characters to a 
variable since the IBS 360/67 computer uses four bytes per 
standard integer word. Names are stored column-wise in 
vectors of three elements. The parameter NH must be set to 
the number of elements required for storage of a name on the 
new system. LBYTE must be set to the number of characters 
per integer word. NH should be set such that NH*L3YTE at 
least 6. Both of the parameters are set in the BLOCK DATA 
subprogram. Data statements initialize operation names in 
the BLOCK DATA subroutine and in all three of the GROOP 
subroutines. All these data statements must be adjusted so 
that they initialize LBYTE characters per element. Format 
statements for write commands with output in format A4 
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should 

should 

should 



be changed. For the most part the output format wAd 
be changed so that w = NH and d = LBTTE. Other uses 
be apparent to the programmer. 



E. POOPPA B X .A Lia miMX. 



There are three potential areas which may cause 

boundary alignment problems. NDP is set in the main program 
to the number of integer words needed to store a real 
variable. NDIR is the number of integer words used for a 
matrix directory. The first possible problem area is the 

relation between NDIR and NDP. If NDIR is not evenly 

divisable by NDP, add the remainder to the equation for NDIR 
in the BLD3K DATA Subroutine. The second possible problem 
area is in the Subroutine LIST. This subroutine contains a 
test to ensure the number of elements reserved for a single 

precision array is evenly divisible by NDP. If NDP is more 

than 2, tie test is not valid. A better test would be 
IF (NP.EQ.1) 1(1) = NR * NC ♦ HOD (NR* NC, ND P) . 

Finally, if NDP is greater than 2, then blank common must be 
adjusted so that there are NDP integer variables preceeding 
the L array in the main program and all subroutines where 
blank common appears. 
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V . 



CONCLUSIONS AND RECOMMENDATIONS 



CAL provides a flexible tool for teaching modern 
sturctural analysis. The general matrix operations make it 
useful in the area of linear algebra as veil as structural 
analysis. The use of CAL is highly encouraged. 

A. RECOMMENDED MODIP IC AT 10 NS 

A graphical display capability should be incorporated in 
CAL as soon as the equipment -becomes available. It would be 
espically helpful to have a display of the element forces 
from the FORCE operation. Additionally, incorporation of 
heat transfer options into the program should be considered. 
The heat transfer capability would be a valuable tool for 
students in the department. 

B. QT-HER APPLICATIONS 

The organization of CAL provides an excellent outline 
for development of other major analytical programs. The 
storage scheme can be simply modified to permit matrix 
storage using the skyline algorithm. Large problems can 
then be accomodated by incorporating existing algorithms for 
an out-of-core equation solver. 
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APPENDIX A 



PROGRAM LISTING 
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APPENDIX B 



USER'S HA NO A L 



This appendix provides details on use of CAL with the 
IBS 360/67 computer at NPS. Section I provides details on 
the command structure. Section II is a summary of commands 
available. Section III provides the job control language 
for executing the program in both the batch and interactive 
modes at NPS. Section IV contains detailed specifications 
for each available command. Finally, section V gives 
direction for solving larger problems with CAL. The 
majority of this appendix, was originally published as 
reference (3) . The author wishes to express appreciation to 
Professor Wilson for permission to use this material. 
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I 



FORM yjD RESTRICTION OF THE LAN3QAGE 



CAL is an interpretive language which is designed to 
manipulate arrays and matrices and to perform several 
standard structural analysis operations. A CAL program run 
involves the reading of the input deck once and executing 
the commands designated by the operation cards as they are 
encountered. Looping operations allow a sequence of 

commands to be executed more than once. 

The input deck is composed of operation cards and data 
cards. The data cards directly follow each operation card 
which requires data (see LOOP operation for exception to 
this). Che operation card -contains the name of the 
operation to be executed, names of arrays associated with 
the opention and integer constants. Examples of the 

general form of this card are 

0P,M1,M2,M3,M4,M5,N1,N2,N3,N4 COMMENTS 
OP, M1,N1,N2 
OP, N1 
OP 

in which OP is the name of the operation to be executed. Mi 
is the name of an array and Ni is an integer. The names of 
OP or Mi are one to eight alphabetic or numeric characters 
to be selected by the user. The first character of a name 
must be alphabetic. The sequence of terms OP, Mi and Ni must 
be separated by commas. Characters following a blank will 
be printed as comments in the output from the program run. 

If an operation attempts to load or generate an array 
which previously existed the program will delete the array 
before the execution of the operation. A new array need not 
be the same size of the old array which had the same name. 
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II. SUMMARY OF COMMANDS 
A. GENERAL COMMANDS 

* indicates a significant change or addition in CAL-NPS 



STAR! 


- 


Initialize for the next problem 




STOP 


- 


Normal termination 




NO 


- 


Temporary supprossion of output 




YES 


- 


lestores output 




LABEL 


- 


Print comments 




READ 


- 


Change logical device for input 


* 


WRITE 


- 


Change logical device for output 


★ 


TIME 


- 


Suppress time printout 


* 


SAVE 


- 


Interrupt a problem 


* 


R ES UM E 


- 


Continue an interrupted problem 


* 


LIST 


- 


List arrays and storage used 


* 


B. GENERAL MA7RTY COMMANDS 




LOAD 


- 


Load user defined matrix 




ZERO 


- 


Create null or unit matrix 




PRINT 


- 


Matrix print operation 


* 


DDP 


- 


Matrix duplication 




ADD 


- 


Matrix addition 




SOB 


- 


Matrix subtraction 




MULT 


- 


Matrix multiplication 




TRAN 


- 


Matrix transpose 




SCALE 


- 


Multiply a matrix by a scalar 




SOLVE 


- 


Solution of linear equations 


* 


DOPSM 


- 


Form sub-matrix from large matrix 




STOSK 


- 


Store sub-matrix in large matrix 




DOPDG 


- 


Form row matrix from diagonal 
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STODS 

MAX 

NORM 

INVEL 

S 33 EL 

LOG 

PROD 

DELETE 



Store row on diagonal 
Evaluate row maximums 
Evaluate matrix norms 
Inverte each term in matrix 
Square root of each term in matrix 
Natural log of each term in matrix 
Evaluate product of all terms in matrix 
Delete matrix from storage 



C. STATIC ANALYSIS OPERATIONS 



NODES 

BOD ND 

BEAM 

TRUSS 

PLANE 

SLOPE 

FRA ME 

LOADI 

LOADS 

ADDSF 

ADDK 

MEMFRC 

D IS PL 

FORCE 



Input joint geometry 

Specify boundary conditions 

Form 3-D beam stiffness matrix 

form 3-D truss sti ffness , matrix 

Form 3 to 8 node plane sti.ffness matrix 

form stiffness matrix from slope/deflection eg. 

?orm 2-D frame stiffness matrix 

Load integer array (connectivity matrix) 

Form load vector 

Form global stiffness and mass matrices 

Add element matrix to global matrix 

Calculate element forces from joint displacements 

Pcint joint displacements 

Evaluate and print member forces 



D. DYNAMIC A NALf SI S OPERATIONS 



FUNG 

STEP 

EIGEN 

DYNAM 

PLOT 



- Generate equal interval time function 

- Integrate dynamic equilibrium equations 

- Evaluate mode shapes and frequencies 

- Evaluate uncoupled equations of motion by 
mode superposition method 

- Line printer plot of joint time history 
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E. LOOP OPERATIONS 



LOOP 


- Start of loop 




NEXT 


- End of loop 




SKIP 


- Conditional skip of 


operations within loop 


F. 


NAMES AVAILABLE FOR USER 


.■SOPB.tmi.UiES- 



0SEH& 

USERB 
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III. J CL FOR EXECUTION 

The i bb control language (JCL) necessary to initiate 
program execution for the OS/MVT (batch) system is: 

//(standard green job card) 

// EXEC CM. 

//SYSIN DD * 



CAL COMMANDS AND DATA CARDS 



/* 

For the CP/CMS (t ime- sharing) system the procedure is: 



use standard login procedure, then 
CP LINK 0099P 191 199 PASS= ESAN 
LOGIN 199 A, P 
CAL 

CAL COMMANDS AND DATA CARDS 

The following information is provided for the user who 
wishes to use his own subroutine or solve larger problems 
then allowed with the default core allocation. The load 
module library for the 0S/M7T operating system is in 
partitioned data set F0099.CAL on DISK02. The load module 
is assigned member name MODULE. The file used as a scratch 
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pad is FTD1F001. File FT02F001 is used as the scratch file 
for the SAVE and RESUME operations. The FT01F001 file is 
currently pointed at the SYSDA disk and FT02F001 is pointed 
at dummy. If for any reason a user desires to use the SAVE 
operation with the 3S/MVT system, both FT01F001 and FT02F001 
must be redifined. See Computer Center TN No. 0141-05, U SER 
LIBRARIES AND SOURCE CODE EDITING UNDER OS for procedures. 
File FT01F001 requires 3404 bytes of storage per element. 
File FT02F001 requires a total of 4 times the dimension of 
the L array plus 8 bytes and should have RECF M=VB S. 

For CP/CMS, the execution routine is 
&TYPE0UT OFF 
CP SET LINELN 130 
GLOBAL T SYSLIB SSPLIB 

FILE3EF 01 DSK-P4 FILE FT&1F001 RECFM VS LRECL 3408 
BLK SIZE 3412 

FILE3EP 02 DSK-P4 FILE FT02F001 RECFM VS LRECL 7290 
BLKSIZE 7294 

FILED EF 04 DSK 

FILEDEF 08 PRT RECFM FBA LRECL 133 (PERM) 

LOAD STORE USER CALH (CLEAR N3MAP XEQ) 

FI01F001 is the scratch pad for the structural analysis. 
FT02F001 is created by the SAVE operation. Both FT01F001 
and FT02FD01 must be available for the RESUME operation. 

Temporary additional storage can be obtained, if needed, 
by typing GETTEMP while in CMS. The GETTEMP execution 
routine is: 

STYPEOUT OFF 
VSET RDYMSG OFF 
CP S2T LINELN 130 
SBEG3 TACK 
VSET BLIP + 
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CP DEFINE T2 314 192 4 
FORMAT T ALL (NOTYPE) 
RELEASE 192 T 
LOGIN 192 P 
LOGIN 191 B,P 
&ENDS I ACK 
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Iff. CAL COMMAND SPECIFICATIONS 



A. GENERAL MATRIX OPERATIONS 



CAL has most of the standard matrix operations plus 
some special array operations which are useful in structural 
analysis. The following is a list of approximately 25 
operations which are used for control and general matrix 
manipulation . 

♦ indicates the formation of a new matrix. A matrix 
previously defined with the same name will be deleted. 

- indicates modification of an existing matrix. 

START 

This operation eliminates all arrays which were perviously 
loaded or generated. 



STOP 

This operation causes normal termination of a CAL program. 



NO 

IES 



These operations are used to selectively suppress output 
from CAL. The NO operations suppresses all printing, except 
diagnostics, until the operation TES is encountered. 
Therefore, in subsequent runs of the same CAL program, 
output which was previously correct need not be reprinted if 
these cards are inserted in the data deck. 



LAB EL, N 1 

This operation will read and print N1 comment cards which 
follow this operation card. Coluun 1 of each card will be 
interpreted as a standard carriage control symbol (i.e. 0 
for double space and 1 for skip to the top of the next 
page) . 
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HEAD, N1 

IRIS OPERATION IS VALID ONLY WITH THE CP/CHS TIME SHARING 
S YS TEM . 

This operation permits the selection of the offline printer 
or the terminal as the input file device. Default is the 
terminal. N1 =4 will read subsequent commands from FILE 
FT04F00 1 on the users p-disk. NT =5 will restore the 
terminal as the input file device. All disk files prepared 
for use with this command should end with either STOP or 
READ, 4. This command will not be executed on the OS/MVT 
(batch) system. 



WRITE, N 1 

THIS OPERATION IS VALID ONLY WITH THE CP/CHS TIME SHARING 
SYSTEM- 

This operation permits the selection of the offline printer 
or the terminal as the output file device. Default is the 
terminal and all error messages will be printed at the 
terminal regardless of the output file device selected. N1 = 
3 selects the offline printer. N1 = 6 restores the terminal 
as the output file device. This command will not be 
executed on the OS/MVT (batch) system. 



TIME 

This operation permits the time printout to be suppressed 
without loss of other output. A second TIME will restore 
the time printout unless the print output is suppressed with 
the NO command. 



SAVE 

The SAVE command creates FILE FT02F001 containing all arrays 
in storage at the time of issuance. Note that saved arrays 
will contain any modification since creation. For instance, 
if a matrix has been reduced by the SOLVE operation, the 
reduced form of the matrix will be stowed in FILE FT02F001. 
Note the requirements discussed previously for use of this 
command with the OS system. 



RESUME 

The RESUME operation reads FILE FT02F001 into memory. Any 
arrays currently in storage will be destroyed. FILE 
FT02F001 must have previously been created on a mass storage 
device using the SaVE operation. Note that FILE FT01F001 
must also be accessible if an interrupted structural 
analysis problem is being resumed. 



LIST 



The LIST operation prints the 
arrays in storage and the amount 



directory information 
of storage used. 



for 
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LOAD, Ml ,N1 , N2 ,N3 

This operation will load an array of real numbers named Ml 
which has N1 rows and N2 columns. The terms of the array 
are punched in row-wise sequence on data cards following 
this operation. If H3 is zero or blank, the cards are 
punched in a format of (8F10.0). If N3 is nonzero, an 
additional card containing the format of the data cards must 
follow this operation ana preceed the data cards. If the 
data is to be 4 numbers per card in field widths of 15, the 
additional cards would contain the following information: 
(4F15.0) . 



+ 

ZERO, H 1 ,N1 , N2,N3, N4 



A real matrix named Ml is created with N1 
columns. The terms in this* matrix will have 
values . 



Ml (1,1) = N3 I 

Ml (I, J) = N4 J 

Therefore, : this operation 
matrices. 



1 y . . • 
— 1 , . . . 

can be 



N 1 
N2 

used 



to form 



rows and N2 
the following 



null or unit 



PRINT, Ml or PRINT, Ml, N1 or PRINT , M 1 , N 1 , N2 or PRINT, Ml, ,N2 

This operation will print the real array named Ml in a 
matrix format of up to 8 columns per line. If N1 is greated 
than zero the operation will read and print N1 comment cards 
wich follow the operation card. N2 is optional but note 
that an extra comma must be used in place of N1 if no 
printed title is desired. The matrix Ml will be printed in 
partitioned form with N2 columns per partition. Lines will 
have N1*13 + 5 characters. N2 defaults to 8, printing 125 
characters per line. The user is cautioned not to overcome 
the capacity of the printing device in use. 



+ 

DU P, Ml, M2 

This operation will form an array named M2 which is 
identical to the array named Ml. 



ADD, Ml , M2 

This operation will replace matrix Ml with the sum of the 
matrices 31 and M2. 



SUB , Ml , M2 

This operation will replace matrix Ml with matrix Ml less 
matrix M2. 
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+ 



MULT, Ml, M2, M3 

This operation generates a new matrix M3 which is the 
product of matrices Ml and M2, or M3 = M1*M2. 



TRAN, Ml, M2 

This operation generates a new matrix M2 which is the 
transpose of matrix Ml. 



SCALE, M 1 , M 2 

This operation replaces each term in the matrix named Ml 
with tne term multiplied by the term M2(1,1) of the matrix 
named m2. 



SOLVE, Ml, M2, N1,N2 or SOLVE, M 1 ,M2, , N2 or SOLVE , M 1 , N1 , N2 or 
SOLVE, Ml, M2, N1 or SOLVE, Ml, M2 or SOLVE, Ml 

If N1 = 3, this operation solves the matrix equation AX=B. 
Ml is the name of the A matrix and M2 is the name of the B 
matrix. Matrix A is t rang ularized- and the results X, are 
stored in M2. 

If N1=1 Matrix A is triangularized only. 

N1=2 For a given B matrix and the a matrix Dreviously 
triangularized, the 3 matrix is replaced by the results X. 

N 1= 3 Matrix A is replaced by its inverse FOR SYMMETRIX 
MATRICES ONLY. 

N2=0 or blank, matrix A is symmetric. If N2 is nonzero 
the matrix A is not symmetric. 

For symmetric matrices, matrix A is factored into the LDL 
form. The diagonal D matrix is stored on the diagonal of A. 
The parameter N2 permits the direct solution of 

non-symmetric systems of equations. If N2 is not equal to 
0, an LO decomposition of matrix A will be performed. No 
direct replacement of Ml by its inverse is availabe for the 
non-symmetric case. Instead, use the ZERO operation to 
create an identity matrix M2 of the same order as Ml. The 
command SOLVE ,M1 ,M2, , N2 will then replace the matrix M2 with 
the inverse of the matrix A. 



DDPSM,M1,M2, N1,N2,N3,N4 

This operation forms a new submatrix named M2 with N3 rows 
and N4 columns from terms within the matrix named Ml. The 



+ 



first tern of matrix M2 
column N2 of matrix Ml, 




will be from row N1 and 
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STOSM,M1,M2,N1,N2 

This operation stores a submatrix named M2 within the matrix 
named HI. The first term of the submatrix M2 will be stored 
at row N1 and column N2 of matrix Ml. The terms within the 
area of Ml in which M2 is stored will be destroyed. 



f 

DUPDG,M1,M2 

This operation forms a new row matrix named M2 from the 
diagonal terms of matrix Ml. 



ST0DG,M1,M2 

This operation stores a row or column matrix named M2 at the 
diagonal Locations of matrix Ml. 



+ 

MAX, Ml, M2 

This operation forms a column matrix named M2 in which each 
row contains the maximum absolute value of the corresponding 
row in matrix Ml. The maximum and its column number is 
printed for each row. 



NORM, Ml, M2, N1 

If N1 = 3, a row matrix named M2 is formed in which each 
column contains the sum of the absolute values of the 
corresponding column of matrix Ml. If N1 /= 0, a row matrix 
named H2 is formed in which each column contains the square 
root of the sum of the sguares of the values of the 
corresponding columns of matrix Ml. 



INVEL, Ml 

This operation replaces each term in the matrix named Ml 
with it’s inverse. 



SQHEL, Ml 

This operation replaces each term in the matrix named Ml 
with the square root of the term. 



LOG, Ml 

This operation replaces each term in the matrix named Ml 
with the natural log of the term. 



Ill 



PROD, Ml , M2 

This operation forms a 1 x 2 array named M2 which contains 
the product of all terms in the matrix Ml. The product, X, 
is stored as two numbers of the form 
X = p*10**E 

in which 32(1) = P and M2 (2) = E, the exponent. 



DELETE, Ml 

This operation will cause the elimination from storage of 
the array named Ml. 
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B. STATIC ANALYSIS OPERATIONS 



The purpose of this series of operations is to form the 
total stiffness and diagonal (lumped) mass matrices for 
systems of two or three-dimensional elements. For 
thr ee-edimensional analysis, there are beam and truss 
elements available. For two-dimensional analysis, there is 
a frame element, a slope/deflection element for beams, and a 
3 to 8 node isoparametric finite element available. 

After the creation of an array containing the 
coordinates of the joints of the system, the specification 
of displacement boundary conditions, the tabulation of 
material and section properties, the mass and stiffness 
matrices are formed for each structural member and placed in 
sequence on low speed storage along with the global equation 
numbers wiich are associated with their stiffness terms. In 
addition, the member force-displacement transformation 
matrices are formed and stored on a separate low speed 
storage cile along with the appropriate displacement 
num bers . 

The NDDES operation is used to specify or generate the 
geometry of the system. The operation BOUND specifies which 
joint displacements exists and assigns internal equation 
numbers to these displacements. Therefore, each joint may 
have from zero to six displacement degrees of freedom. 
Tables of material and section properties for the various 
members are loaded and printed as standard arrays of 
information. 

A special operation, • ADDSF, is used for the direct 
addition of element stiffnesses to form the total stiffness 
and diagonal mass matrix of the system. The ADDK operation 
may be used to add individual elements into the total system 
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matrices. The LOADS operation specifies the concentrated 
joint loads for all load conditions. After the direct 
solution for joint displacements due to static or dynamic 
loads, the member forces can be evaluated useng the FORCE 
operation. Individual member forces can be evaluated using 
the MEMFR3 operation. The DISPL operation is used to print 
the displacements in joint number order. 
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NODES , M 1,1 1 

The cards following this operation provide information for 
the creation of a N1 by 3 array which will contain the 
coordinates for all the joints in the system. Bhere 



Ml = Name of new coordinate array to 


be loaded. 


N1 = Number 


of joints (or nodes) in 


the system. 


The following sequence of cards punched 
format must follow this operation. 


in a ( 21 5, 3F1 0 . 0) 


Columns 


Contain 




1 - 5 


Node number selected by 


user 


6-10 


blank 




11 - 20 


X-coordina te 




21 - 30 


I-coordinate • 




31 - 40 


Z-coordina te 




Nodes cards may 
must be defined 
last definintion 
be terminated wi 


be supplied in auiy order; however, all nodes 
. If nodes are defined more than once the 
will be used. This sequence of data must 
th a blank card. 



z 




POSITIVE COORD I HATE SYSTEMS 
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BOUND, Ml 



This opention specifie 
for the structural sys 
operation. Where 


Ml = Name 
gen erat ed. 


of bo 


This od 
the fol 


eration 
low ing 


is foil 
inf orma ti 


Columns 




Conten 


1 


- 5 




Node n 
series 
ment s 


6 


- 10 




Node n 
series 


1 1 


- 15 




X-tran 


16 


- 20 




Y-tran 


21 


- 25 




Z- tran 


26 


- 30 




X-rota 


31 


- 35 




Y-rota 


36 


- 40 




Z-rota 


41 


- 45 




Node n 
Condit 


A transit tion or rota 
undifined displacemen 

displacements to be eva 



If a node boundary 
displacema nts at that 
supplied in any order 
specified iiore than one 
sequence of data must b 

The selection by the 
displacements requires 
stiffness procedure. D 
have no stiffness asso 
considered to be undi 
develop an equilibrium 
total number of nonzero 
size of the total sti 
ADDSF operation. 



s the displacements which are nonzero 
tern of joints specified by the NODES 



undary condition cade array to be 



owed by a series of cards containing 
on punched in a (815) format. 

ts 

umber for the first node in a 
of nodes with identical displace- 
pecif ication. 

umber for the last node in the 



slation. 



slat ion 

slation 

tion 

tion 

tion 



umber increment used to generate 
ions for additional nodes. 



tion equals: (a) zero for zero or 
ts, or (b) one for nonzero 
luated by other operations. 



condition is not specified, all 
node are assumed zero. Cards may be 
. If node boundary conditions are 
e, the last difinition is used. This 
e terminated by a blank card. 



user of which 
an understand 
isplacements deg 
ciated with the 
fined since it 
equation for 
displacements 
ffness matrix t 



nodes have nonzero 
ing of the direct 
rees of freedom which 
displacement must be 
is not possible to 
that direction. The 
specified will be the 
o be defined by the 
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BEfcM,M1,M2,M3,M4 



This operation calculates the element stiffness, mass and 
force-displacement transformation matrices for 3-D beam 
members. These arrays are stored in sequence on low speed 
storage to be used by other operations where: 

Hi is the name of the beam element group 

M2 is the name of the coordinate array 

M3 is the name of the boundary condition array 

M4 is the name of the arra which contain beam 

properties and has been loaded by the standard 
matrix LDAD operation 



One card for each beam in this group of beam elements must 
follow this operation. The beam cards are punched in (515) 
format, where 



Columns 
1 - 5 

6-10 
11 - 15 

16 - 20 
21 - 25 



Contain 

Beam identification number 
Node number‘1 
Node number J 
Node number K‘ 

Beam property number NP 



This sequence of cards must be terminated with a blank card. 

The material and geometric properties for each element are 
give in the M4 array in the following order: 



a 4 

M4 
M4 
H4 
M 4 
M4 
H4 



(N ? , 1 
(NP, 2 
(NP, 3 
(NP, 4 

<;»?,§ 
(NP , 6 

i NP, 7 



Axial area of member, A 

Torsional Moment of Inertia, J 

Moment of Inertia about axis 2, I 

Moment of Inertia about axis 3 , I 

Modulus of Elasticity, E 

Shear Modulus, G 

Mass per unit length of beam 



where NP is the specific material property number specified 
in columns 21 - 25 of the beam card. The local sign 
convention is given in the following figure. 
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2 





local beam reference system. 
Axes 1 and 2 are in the 
plane defined by plane 
i, J and k nodes. Axis 1 
is defined by line i-J . 

Axis 3 is perpendicular 
the i , J , k plane. 




DEFINITION OF POSITIVE BEAM FORCES 
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TRUSS, M 1,.12, M3,M4 



This operation forms the element stiffness, mass and 
force-displacement transformation matrices for 3-D truss 
■embers. The arrays are stored on low speed storage in 
sequence and will be used by other structural operations. 



Ml is the name of this groups of truss members 

M2 is the name of the coordinate array 

M3 is the name of the boundary condition array 

M 4 is an NP by 3 array of section properties in 

which NP is the number of different section 
properties and 

M4 (N ? , 1 ) = the cross-section area, A 

M4 [NP, 2) = the Modulus of Elasticity, E 

M4/NP,3f = the mass per unit length of the member 

This matrix can be loaded by the matrix load operation. 

This operation is followed by one card per truss member in 
(415) romat with the following information: 



Columns 
1 - 5 

6-10 
11 - 15 

16 - 20 



Content 

Truss member identification number 

Joint number I 

Joint number J 

Section property number, NP 



This operation must be terminated by a blank, card. 



+ 4 * 

LDADI,M 1,N1,N2,N3,N4 or LO ADI , M 1 , N 1 , N2 , ,N4 or 
LOADI,M1,N1,N2,N3 or LO ADI , Ml , N 1 , N2 

This operation will load an inteaer array named Ml which has 
N1 rows aid N2 columns. The terms of the array are punched 
in row-wise sequence on data cards which follow this 

operation. If N3 is zero or blank, the data mu?t be punched 
in (16I5| format. If N3 is nonzero, an additional card 
containing the format of the data cards must follow this 

operation and preceed the data cards. N4 is optional but 

note that an extra comma must be used in place or N3 if the 
user does not supply a format. The matrix Ml will be 

E rinted in partitioned form with N4 columns per partition. 

ines hava (N4 + 1) *5 characters. N4 devaults to 20 causing 
125 characters to print per line. The user is cautioned not 
to overcome the capacity of the printing device in use. 
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PLANE,M1,1 2 , M3,M4 ,N 1, N2 

This operation calculates the element stiffness, mass and 
str ess- displacement transformation matrices for 3 to 8 nod 
isoparametric elements. (Y-Z plane only) . These arrays ar 
stored in sequence as a group on low speed storage to be 
used later by other operations (i.e., ADDS? and FORCE). The 
arguments are defined as 

Ml is the user defined name of the element group 

M2 is the name of the joint coordinate array 

M3 is the name of the boundary condition array 

M4 is the name of the array which contains the 

material properties of the elements (one row per 
different material) where 



MU (N? , 1 
M4 (NP, 2 
M4 (NP, 3 
M4 (N? , 4 



= Modulus of Elasticity, E 
= Poissons Ratio, v 
= Thickness of the element 
= Mass density of the element 

NP is the material identification 



number 



and 



N1 and N2 are the number of 
s directions respectively. 



integration points in the r 



One card for each 3 to 8 node element in the group must 
follow the operation card. The card is punched in a 
( 13 15, 6F5. 0) format and contains the following information: 

Contain 

element identification number 
Node number N1 
Node numoer N2 
Node number N3 
Node number N4 
Node number N5 
Node number N6 
Node number N7 
Node number N8 

Material identification numb 
Natural coordinate of stress 
Natural coordinate of stress 
Natural coordinate of stress 
Natural coordinate of stress 
Natural coordinate of stress 
Natural coordinate of stress 



Columns 


1 


— 


5 


6 


- 


10 


1 1 


- 


15 


16 


- 


20 


21 


- 


25 


26 


— 


30 


31 


- 


35 


36 


- 


40 


4 1 


- 


4 5 


46 


- 


50 


51 


- 


55 


56 


- 


50 


61 


- 


65 


66 


- 


70 


71 


- 


75 


7 6 


— 


30 


N4 thro 


ugi 


. N 



er, NP 
output rl 
output si 
output r2 
output s2 
output r3 
output s3 



are optional. The aidside nodes, if present, 
must be'#ithin center half of side. The local numbering 
system for the element is shown in the following figure. 
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I SOP ARAM ETHIC ELEMENT 

Stresses will be printed by the FORCE operal 
points defined in columns 51 throuah 80. 
defined as follows: 



*1 




■ T (Dl 

l xx 


I CM ( 
Pm Cm 




X U) 
> 


3 

F. 




^(2) 


4 

F- 


— 


t ^2) 


5 

F 6 

F_ 




t ^2) 

T Xy (3) 


7 

F 8 

F. 


• 


t ^3) 




"(3) 


9 







ion at the three 
The forces are 
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SLOPE, Ml 



rhis operation forms a 4 x 4 stiffness matrix. Hi for a beam 
or column member from the classical slope-deflection 
equations. The properties of the member are defined on one 
immediately following the operation. This second card 
1 (3F10.0) format and contains the following 



car d 

is punched in 
information : 
Columns 

1 - 10 
11 - 20 
21 - 30 



Contain 

Moment of Inertia, I 
Modulus of Elasticity, 
Length of Member, L 



The sign convention is defined as follows: 




The memner forces are defined in terms of joint 
displacements by the following slope deflection equations. 



M 2 



ML 

L 

ML 

L 



40, + 20 2 - 
20, + 40 2 - 



Yi--V 3 ~ 



a/,+a / 2 



6 _ 

L 

6 _ 

L 



('■ " "=)] 
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or in matrix form 







4 2 -A A 

L L 




m 7 


El 


2 4 -A 1 

L L 


*2 


V\ 


L 


6 6 12 12 


v \ 






L L L 2 L 1 

6 6 12 12 

L L L 1 L 1 


v 2. 



or symbolically F = KA Where K is the 4 x 4 stiffness 
matrix defined by the name HI. 



+ * 

FRAME, HI, M2 

This operation forms the 6x6 stiffness matrix Ml for the 
two-dimensional frame member shown below. 
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£he properties of the member are defined on one card 
immediately following the FRAME operation card. This second 
card is punched in a (8F10.0) format and contains the 
following information: 



Coluans Contain 

1-10 Axial area, A 

11 - 20 Modulus of Elasticity, E 

21 - 30 Moment of Inertia, I 

31-40 Xi 

41-50 Yi 

51 - 60 Xi 

61 - 70 ij 




M2 is a 3 x 6 
is based on 
shown below. 



force-displacement transformation matrix which 
the positive difinition of the element forces 
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These forces can be calculated from the following matrix 
equation, with the MEMFRC operation. 



M, 

Mj 

P 



M2 



U, 

Uy. 

Uu. 

Uyj 

U.J 



4 

LOADS, Ml, M2, N1 



This operation forms a load matrix named Ml of N1 columns 
(N1 load conditions) where M2 is the name of the boundary 
condition array generated by the operation BOUND. This 
operation is followed by a series of cards - one for each 
loaded joint for each load condition. These cards are 
punched in (2I5,6F10.0) format as follows: 



Colua ns 
1 - 5 

6-10 
11 - 20 
21 - 30 
31 - 40 
41 - 50 
51 - 60 
61 - 70 



Contain 
Joint number 
Load condition number 
Load in X-direetion 
Load in X-direction 
Load in Z-direction 
Moment about X-axis 
Moment about X-axis 
Moment about Z-axis 



This series of cards must be terminated by a blank card. 



♦ ♦ ♦ 

ADDSF,M 1 or ADDSF,M 1 , M2 

This operation forms the total stiffness matrix named Ml and 
a lumped lass matrix named M2 for the structural system from 
the element stiffness and mass matrices which are stored on 
low speed storaae. These matrices can be printed with the 
PRINT operation. If M2 is not specified, the row mass 
matrix M2 will not be formed. 



ADDK,M1,M2,H3,N1 

This operation adds the element stiffness matrix named M2 to 
the total stiffness matrix named Ml, where Ml was previously 
defined and initially set to zero. M3 is the name of the 
integer array in which the column number HI contains the row 
or column numbers in the total stiffness matrix where the 
element stiffness terms are to be added. 
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ME1FRC,M1,M2,M3,M4,N1 

This operation multiplies the element stiffness matrix named 
Ml by the joint displacement matrix named M2. M3 is the 
name of the integer array in which the column number N1 
contains the row numbers in the displacement matrix, M2, 
which are to be multiplied by the element stiffness (or 
force-displacement) matrix., HI. The results of this 
multplication are stored in the array named M4. 



DIS PL, M 1 , M2 

This operation prints the displacement named Ml in joint 
sequence order, where M2 is the name of the boundary 
condition array. 



FORCE, Ml, M2, M3 or FOR CE, M 1 ^ M2 

This operation calculates the - member forces for a group of 
elements in whcih 

Ml is the name of the element group 
M2 is the displacement matrix 

M3 is the name of the mat-rir in which the forces are 
stored in the order calculated. If this array is not 

specified, the element forces will be printed only and will 
not be retained in storage. For the TRUSS element, only the 
member axial force, F, will be calculated for each member. 
For the 3EAM element, eight forces will be printed with 
reference to the positive definiton shown in the BEAM 
operation. 
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c. 



DYNAMIC ANALYSIS OPERATIONS 



The following operations were designed to evaluate the 
dynamic response of structures subjected to arbitrary 
time-dependent loads. If these operations are used in 
connection with the standard matrix operations and the 
sturctural analysis operation, a dynamic analysis is a 
relatively simple procedure. The user has the option of 
using the mode superposition method or a direct step-by-step 
integration of the dynamic equations of motion. The user 
may examine the spectra of both input loading and calculated 
displacements. In additiorf, the contributions of the 
individual modes may be evaluated and compared. 

The most common and convenient form for time-dependent 
data to be specified is as straight line segments between 
given time points. Therefore, an operation which generates 
values at equal intervals is necessary. Another common 
characteristic of time-varying loads on structures is that 
it is normally possible to represent the loads at all points 
on the structure by the product of two matrices, a column 
matrix indicating the spacial distribution of loads times a 
row matrix which indicates the values as a function of 
various times. If a more complicated loading is required, 
it is possible to perform more analyses, each within the 
restrictions of the program, and then add the results of 
each analysis. 

The following operations have been added for the major 
purpose of performing dynamic analysis. 
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FUNG, Ml, M2, M3, N1, N2 ' * 

This operation generates a matrix named H2 which contains 
values at equal intervals of the function specified in the 
array named Ml. The array Ml must be a 2 by k array of the 
form: 



M 1 



h *2 

*1 <2 



which numerically 
below: 



*3 \ 

S \ 

represents a function of the form shown 




The time interval t is specified in the 1x1 matrix named 
M3. N 1 specifies the total number of values to be 
generated, and is the number of columns in M2. If N2 = 0. 
the array M2 will be a 1 x N1 row matrix in which the first 
value will be f . If N2 0, the array M2 will be a 2 x N 1 
matrix of the following form: 



t x t x + at t + 2At ... 

f(t 1 +At) f(t x + 2At ) 



M2 * 



+ 

STEP, Ml , M2 , M3 ,M4 , M5,M6, M7,M8, N1 , N2 



This operation calculates the dynamic response of 
sturcturaL system usnig direct step-by-step integration o 
-the following linear matrix equation or motion: 

HO. + CO ♦ KU = R (t) = PF (t) 

Where 

Ml is the name of the N x N stiffness matrix K 

M2 is the name of the H x N mass matrix M 

M3 is the name of the N x N damping matrix C 

M 4 is the name of the N x 3 initial condition matrix 

0 in which: 



0 (1,1) is a vector of displacements 0 

0 (1,2) is a vector of velocities 0 

U (1.3) is a vector of accelerations U 

M5 is the name of the N x N2 matrix of calculated 

displacements in which column "i" represents the 
displacements at time i*N1*At 
M6 is the name of the N x • 1 load distribution matrix 
P 

M7 is the name of the 1 x k row matrix representing 
the load multipliers at equal time increments 
F, where k = N2/N1 

M8 is the name of the 1x1 matrix containing ^t 
N 1 is the output interval for the displacements 
N2 is the total number of displacement vectors to 
be calculated. 



The total time for which results will be calculated by this 
operation is N1*N2*At. This operation must be followed with 
one data card in (3F10.0) format containing the following 
information: 

Columns Contain 

1 - 10 DELTA 

11-20 ALPHA 

21 - 30 THETA 

Different values of delta, alpha and theta will allow the 
user to select different methods of step-by-step 
integration. The following table lists some possibilities: 









DELTA 


ALPHA 


THETA 


Newmar ks 


Average Acce 


leration 


1/2 


1/4 


1.0 


Linear Ac 


celeration 




1/2 


1/6 


1.0 


Wilson's 


Theta Metnod 


(low damping) 


1/2 


1/6 


1.42 


Wilson 's 


Theta Method 


(high damping) 


1/2 


1/6 


2.0 
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EI3EN,M 1,M2,M3,N1 



This operation solves the following eigenvalue problem: 

= M#A 

In which the N x N, symmetric, positive-semidef inite matrix 
K is named Ml. The matrix M is a diagonal matrix of 

nonzero, positive- terms designated by M3. The matrix M3 

must be a row or column matrix containing only the diagonal 
terms of M. The eigenvalues, A , are stored m matrix M3. 
The eigenvalues are ordered m numerically increasing order 
and the eigenvectors,# , are stored in the corresponding 
columns of the matrix M2. The number N1 specifies the 
approximate number of significant figures of the 

eigenvalues. If N1 is zero or blank, 4 figure accuracy will 

be used. The maximum accuracy possible is 16 figures. The 
use of more than 12 figure accuracy is not recommended. 

The program reduces the problem to standard eigenvalue form 
by the following transformation 

where 



= m t 


K 


m 


t 


M 




= m 


m 



in which 






1 / 




The calculated mode shapes, $ 



are normalized as follows: 




p K, <j> A 



rhe program uses 
to solve for all 



the standard Jacodi d iagonaliza tion method 
eigen values and eigenvectors. 
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♦ 

DIN AM,M 1,1 2, M3,M4 ,M5, MS ,N1 



This operation evaluates the followi 
second order differential equations ass 
superposition method for the dyna 
structural system. 

X; + 2 XiWLXi + U £ X; . * I? ft) 

Ml is the name of a row or column matr 
N terms (frequencies in rad/sec) . 
row or column matrix which contain the 
modal damping to critical damping) . 

The generalized time-varying forces P: 
directly but are evaluated from 
information. The forces for all mod 
specific times by the program from 
equation: 

P = p*f = M3* M4 



In 
1 x 
arc 
arr 
nec 



which £ is a specified N x 
N 1 row matrix which will 
ay named M4. The array M4 
ay described under the 
essary to use FUN3 before 



1 vector 
be aener 
is the s 
operation 
the DINAH 



M5 is the name of the N x N1 arra 
generalized displacement X^ (t) . 

M6 is the name of the 1 x 1 array whi 
increment associated with the generaliz 



N1 is the number of displacements to be 



The method of integration used is exa 
secments. 



ng set of uncouples 
ociated with the mode 
mic analysis of a 



i = 1 to N nodes 

ix which contains the 
M2 is the name of a 
N A; terms (ratio of 



(t) 


are 


not sp 


ecif ied 




more 


fund 


amental 


es 


are 


evalua 


ted at 


the 


following 


mat rix 



named M3, and f is a 
ated from the 7 by k 
ame form as the input 
F0N3. It is not 
operation. 

y which contain the 



ch contains the time 
ed displacements. 

generated. 

ct for straight line 
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PLOT, Ml ,N1 

This operation will prepare a printer plot of selective rows 
of the matrix named Ml. N1 is the number of rows of HI 
which will be plotted by this operation. This operation is 
foLlowed by N1 cards in (1A1,I4) format with the following 
information: 

Columns Contain 

1 Plot symbol - any keypunch symbol 

2 - 5 Row number to be ploted 

The program automatically searches the information to be 
plotted for the maximum and minimum values. The difference 
in these numbers divided by 120 spaces is selected as the 
plot scale. 
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D. LOOPING OPERATIONS 



CAL has a five level looping ability. The first 

operation is LOOP and the last operation is NEXT. 

Operations withing CAL are normally executed as they are 
encountered. If the operation requires data, the data cards 
follow the operation card. In the case of looping, however, 
alL operation cards from the first LOOP card to the last 
NEXT card are stored within the computer before they are 
executed; therefore, if operations within the loops regure 
data, the data cards must be, supplied in the order required 
after the last NEXT operation.. If an error is encountered 
while executing in a loop, the entire matrix of loop 
commands is deleted and the user is given the opportunity to 
try again. Matrices that have been modified by operations 
successfully completed while in the loop remain modified. 
After all loops are executed the computer storage required 
for these operations is automatically released by the 
program. The looping operations are; 
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LOOP, N 1 



N1 is the number of times the loop is to be executed. 
Associated with each LOOP operation there must be a 
corresponding NEXT operation which signifies the end of the 
Loop and the return of the control to the beginning of the 
loop. Tie following is a possible series of looping 
operations . 



LOOP, 5 <- 



LOOP, 2* 



NEXT- 



Second level loop 
executed twice 
(total of 10 times) 



LOOP, 4* 



NEXT«— 

NEXT--*- 



Second level «1 o o p 
executed 4 times 
(total 20 times) 



First level loop 
executed 5 times 



data for all operations within all loops 



NEXT, Ml or NEXT 

The operation NEXT signifies the end of a loop. It is 
apparent *hich LOOP and NEXT cards are associated if there 
are an egual number of each. The operation NEXT, Ml will 
cause the loop to terminate if the first term in the matrix 
naaed Ml is negative. 



SKIP, Ml ,N1 

This operation will cause the skip of the next N1 operations 
if the first term in the matrix named Ml is negative. This 
level of looping. 



E. OSER DEFINED OPERATIONS 



OSERA and OSERB 



These names 
programmed 
operations 
organization 



are reserv 
by the j 
is is nec 
of CAL. C 



ed for 
ser. 
essary 
hapt er 



operations to be defined and 
In order to program these 
to understand the internal 
III contains details. 
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V 



LARGE PROBLEM S 



CAL is designed as an educational tool. It does not 
take advantage of banding and symmetry in matrix storage. 
Larger problems can be solved by increasing the dimension of 
the L array but a general purpose program that makes maximum 
advantage of out-of-core storage and takes advantage of 
banding and symmetry for in-core matrix storage is probably 
a better choice. With the above disclaimer, to increase 
problem size capability, increase the dimension of the L 
array and change the value of MAX to the new dimension size 
in the follosing: 

C MAIN PROGRAM 

C SET PROGRAM CAPACITY 

COMMON NT0T,NDP,L (6000) 

MT 0T=5 000 
NDP=2 

CALL SETIME 
CALL CALI 
STOP 
END 

SETIME is an OS routine to initialize the CPU timer and 
should not be used for CP/CMS programs. With the dimension 
of the L array as above, the program currently executes in 
144k bytes for OS and 256k bytes for CP/CMS. The region 
necessary for execution will increase about eight times the 
increase in the L array. 
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